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Abstract. During the several months of the outburst of 
Mkn 501 in 1997 the source has been monitored in TeV 
7-rays with the HEGRA stereoscopic system of imag- 
ing atmospheric Cherenkov telescopes. Quite remarkably 
it turned out that the shapes of the daily 7-ray energy 
spectra remained essentially stable throughout the en- 
tire state of high activity despite dramatic flux variations 
during this period. The derivation of a long term time- 
averaged energy spectrum, based on more than 38,000 
detected TeV photons, is therefore physically meaning- 
ful. The unprecedented 7-ray statistics combined with 
the 20% energy resolution of the instrument resulted in 
the first detection of 7-rays from an extragalactic source 
well beyond 10 TeV, and the first high accuracy mea- 
surement of an exponential cutoff in the energy region 
above 5 TeV deeply into the exponential regime. From 
500 GeV to 24 TeV the differential photon spectrum 
is well approximated by a power-law with an exponen- 
tial cutoff: dN/dE = N (E/l TeV)~ Q exp (-E/E ), with 



N = (10.8 ±0.2 stat ± 2.1 sya ) • 10" u ena- 
ct! = 1.92 ± 0.03 sta t ± 0.20 sys , and E = (6.2 ± 
0.4 sta t (—1-5 + 2.9) sys ) TeV. We summarize the methods 
for the evaluation of the energy spectrum in a broad dy- 
namical range which covers almost two energy decades, 
and study in detail the principal sources of systematic er- 
rors. We also discuss several important astrophysical im- 
plications of the observed result concerning the production 
and absorption mechanisms of 7-rays in the emitting jet 
and the modifications of the initial spectrum of TeV radi- 
ation due to its interaction with the diffuse extragalactic 
background radiation. 
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1. Introduction 
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Mkn 501, an active galactic nucleus (AGN) at a redshift 
z = 0.034, was discovered several years ago as a faint 
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ity, unique in both its strength and duration. The TeV 
emission of the source from March to September 1997 was 
characterized by a strongly variable flux. It was on aver- 
age more than three times larger than the flux of the Crab 
Nebula, the strongest known persistent TeV source in the 
sky. Fortunately, the time period of the outburst coincided 
with the source visibility windows of several ground-based 
imaging atmospheric Cherenkov telescopes (IACTs) de- 
signed for the detection of very high energy (VHE) cosmic 
7-rays. Thus almost continuous monitoring of Mkn 501 in 
TeV 7-rays with several IACTs (CAT, HEGRA, TAC- 
TIC, Telescope Array, Whipple) located in th e Nor thern 
Hemisphere was possible (e.g. Protheroe et al. 1997). 



The observations of Mkn 501 by the HEGRA stereo- 
scopic IACT system during this long outburst made a de- 
tailed study of the temporal and spectral characteristics 
of the source possible, based on an unprecedented statis- 
tics of more than 38,000 TeV photons (Aharonian et al. 
1999a; hereafter Paper 1). The "background-free" detec- 
tion of 7-rays , with an average rate of several hundred 
7-rays per hour (against ~ 20 background events caused 
by charged cosmic rays), allowed us to determine statis- 
tically significant signals for < 5 minute intervals during 
much of the 110 h observation time, spread over 6 months. 
Moreover, it was possible to monitor the energy spectrum 
of the source on a daily basis. Within the errors the en- 
ergy spectrum maintained a constant form over the range 
from 1 TeV to 10 TeV. This was the case even though the 
flux varied strongly on time scales < 1 day. We believe 
that this is an important result, and it was to some extent 
unexpected. 

The diurnal spectra exhibit a power-law shape at low 
energies (between 1 TeV and several TeV), with a gradual 
steepening towards higher energies (Paper 1). Such a spec- 
tral form could not be unequivocally ensured in the first 
analysis which was performed during the period of activity 
of the source, since the systematic errors of the recently 
commissioned stereoscopic system of HEGRA were not 
well studied at this time. As a consequence, the energy 
spectrum could not be determined more precisely than 
implied by a power law fit (Aharonian et al. 1997a), even 
though the tendency for a gradual steepening of the ob- 
served spectra was noticed (Aharonian et al. 1997b). The 
results of Paper 1 and the new results in a broader energy 
interval presented below are based on detailed systematic 
studies (see Paper 1 and Konopelko et al. 1999| ), and ex- 
tend and supersede these previous results. This allowed us 
to come to the definite conclusion that the spectrum de- 
termined in the energy region from 1 to 10 TeV steepens 
significantly (Paper 1). A similar ten dency has been found 
also by the Whipple (Samuelson et al. 1998; ), Telescope Ar- 
ray (Kajino et al. 1999, private communication), and CAT 
(Djannati-Atai et al. h999|) groups. Independent spectral 



Apart from its astrophysical significance, the con- 
stancy of the spectral shape has the important practical 
consequence that it allows to measure the spectrum with 
small statistical errors also in the energy regions below 1 
TeV and above 10 TeV. Indeed, the low photon statistics 
of the detector in both "extreme" energy bands (towards 
low energies basically due to the decrease of the detector's 
collection area; towards high energies due to the steep pho- 
ton spectrum) can be drastically increased by using the 
data accumulated over the whole period of observations. 

In Sect. 2 we describe the HEGRA stereoscopic system 
and the specific form of the data analysis, based on Monte 
Carlo simulations of both, the air showers and the detec- 
tion system. The data sample is the same as in Paper 1 
and is described in Sect. 3. A detailed study of the sys- 
tematic errors in the spectrum derivation is contained in 
Sect. 4; most of this methodology was actually developed 
in the context of the Mkn 501 data analysis. We believe 
that this is the first study of this kind. The experimental 
results are then presented in Sect. 5, whereas Sect. 6 at- 
tempts a first discussion. This discussion concentrates on 
the 7-ray results of Sect. 5 and what one can learn from 
them alone. A multiwavelength analysis is clearly outside 
the scope of this paper. Our conclusions are contained in 
Sect. 7. 

Readers only interested in the astrophysical results, 
should skip Sect. 2-4 and proceed to Sect. 5. 



2. The HEGRA system of imaging atmospheric 
Cherenkov telescopes 

2.1. The HEGRA Cherenkov telescope system 

The VHE 7-ray observatory of the HEGRA collaboration 
consists of six imaging atmospheric Cherenkov telescopes 
(IACTs) located on the Roque de los Muchachos on the 
Canary island of La Palma, at 2200 m above sea level. 
A prototype telescope (CT1) started operation in 1992 
and has undergone significant hardware upgrades since 
then. This telescope continues to operate as an indepen- 
dent instrument. The stereoscopic system of Cherenkov 
telescopes consists of five telescopes (CT2 - CT6), and 
has been taking data since 1996, initially with three and 
four telescopes, and since 1998 as a complete five-telescope 
system. Four of the telescopes (CT2, CT4, CT5, CT6) are 
arranged in the corners of a square with roughly 100 m 
side length, and one telescope (CT3) is located in the cen- 
ter of the square. During 1997, when the data discussed in 
this paper were taken, CT2 was still used as stand alone 
detector. 

The telescopes have an 8.5 m 2 tessellated reflector, fo- 
cusing the Cherenkov light onto a camera with 271 pho- 
tomultipliers (PMTs), covering a field of view of 4.3° in 
diameter. A telescope is triggered when the signal in at 
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order to trigger the CT system and to initiate the readout 
of data, at least two telescopes have to trigger simulta- 
neously. Typical trigger rates are in the 10-16 Hz range. 
The PMT signals are digitized and recorded by 120 MHz 
Flash- ADCs . The telescope hardware is described by Her- 
mann (1995); the trigger system and its performance are 
reviewed by Bulian et al. ( 1998| ). 



2.2. Reconstruction of air showers with the HEGRA 
I ACT system 

The routine data analysis (see Paper 1 for details) in- 
cludes a screening of data to exclude data sets taken at 
poor weather conditions or with hardware problems. In 
particular, the mean system trigger rate proved to be a 
sensitive diagnostic tool. Reconstruction of data involves 
the deconvolution of Flash- ADC data (Hess et al. 1998| ), 
the calibration and fiat-helding of the cameras, the deter- 
mination of Hillas image parameters, and the reconstruc- 
tion of geometrical shower parameters based on the stereo- 
scopic views of the air sh ower obtained with the differen t 
telescopes (Daum et al. 1997 , Aharonian et al. 1997a ). 



The characteristic angular resolution for individual 7-rays 
is 0.1°; by sophisticated procedures 7-ray sources can be 
locat ed with sub-arcminute precision (Piihlhofer et al. 



1997| ) 



The separation of hadronic and electromagnetic show- 
ers is based on the shape of the Cherenkov images, in 
particular using the width parameter. The width of each 
image is normalized to the average width of a 7-ray image 
for a given impact parameter of the shower relative to the 
telescope, and a given image intensity. Here, impact dis- 
tances are obtained from the stereoscopic reconstruction 
of the shower geometry. Cuts are then applied to the mean 
scaled width obtained by averaging the scaled width values 
over telescopes. For a point source both the pointing infor- 
mation and the image shape information are used. Each 
of them provides a cosmic-ray background rejection of up 
to 100. 

The reconstruction of the energy of air showers is based 
on the relation between the shower energy and the im- 
age intensity {size) at a given distance from the shower 
axis (Aharonian et al. 1997c). This relation is tabulated 
based on Monte Carlo simulations, with the zenith angle 
of the shower as an additional parameter. The distance 
between a given telescope and the shower core is known 
from the stereoscopic reconstruction of the shower, with a 
typical precision of 10 m or less, for not too distant show- 
ers. The energy estimates from the different telescopes are 
then averaged, taking into account the slightly different 
sensitivities of the telescopes. These sensitivities are cal- 
ibrated to 1% by comparing the light yield in two tele- 
scopes for events with cores halfway between the two tele- 



olution of 15% to 20%, depending on the selection of the 
event sample. 

2.3. Monte Carlo simulations of air showers and of the 
telescope response 

Any quantitative analysis of IACT data has to rely on 
detailed Monte Carlo simulations to evaluate the detection 
characteristics of the instrument. 

The simulation of air s hower s and of Cherenkov light 
emission (Konopelko et al. 1999 ) includes all relevant el- 
ementary processes. On their trajectory to the detector, 
photons may be lost by ozone absorption, Mie scatter- 
ing, and Rayleigh scattering. Atmospheric density pro- 
files, ozone profiles and aerosol densities have been checked 
against local experimental data where available (e.g. Hem- 
berger [1998^ . 

On the detector side, the simulations include the wave- 
length dependence of the mirror reflectivity, of the light 
collection system, and of the PMT quantum efficiency. 
The point spread function of the mirror system is mod- 
eled after measurements of images of bright stars. The 
readout electronics is simulated in significant detail. PMT 
output waveforms are modeled by superimposing the re- 
sponse to single photoelectrons, with their relative timing 
and amplitude smearing. These signals are then sampled, 
quantitized, and fed into the same analysis path as reg- 
ular Flash-ADC data. The simulation includes the mea- 
sured saturation effects both in the PMT/preamplifier and 
in the Flash- ADC. Details concerning the Monte Carlo 
simulation used here, the performance of the system, and 
the comparison with experimental data are described by 
Konopelko et al. (|1999[). 



3. The Mrk 501 data sample 

The analysis is based on the same data sample as used in 
Paper 1, corresponding to a total exposure time of 110 h, 
between March and October of 1997. Mrk 501 was ob- 
served in the so-called wobble mode, with the source posi- 
tioned ±0.5° in declination away from the optical axis of 
the telescopes, alternating every 20 min. For background 
subtraction, an equivalent region displaced by the same 
amount in the opposite direction is used. The separation 
of 1° of these on-source and off-source regions is large com- 
pared to the angular resolution of the telescope system. In 
total, the sample comprises about 38,000 7-ray events. 

To select 7-ray candidates, the same loose cuts were 
applied as in Paper 1. In particular, the reconstructed 
shower direction had to be within 0.22° from the source, 
and the mean scaled width parameter had to be less than 
1.2. Events were accepted up to a maximum impact pa- 
rameter of 200 m from the central telescope CT3. Events 
with larger impact parameters frequently suffer from trun- 
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shower parameters. Since the analysis presented in the fol- 
lowing emphasizes the control of systematic errors, it was 
felt that in this case the advantage of having clean events 
and a well-defined effective area at large energies - basi- 
cally all events within the 200 m radius trigger above a few 
TeV - outweighs the gain in statistics which could have 
been achieved by accepting all events [] For the follow- 
ing analysis a software threshold of two or more triggered 
telescopes, each with at least 40 recorded photoelectrons 
was used. 

Properties of this Mrk 501 data sample were examined 
in detail in Paper 1. 

4. Determination of energy spectra and sources of 
systematic errors 

Compared to single IACTs, stereoscopic IACT systems 
permit drastic reduction of systematic errors, in partic- 
ular for tasks like the precision determination of energy 
spectra. Using the redundant information provided by the 
multiple views, essentially all relevant characteristics, such 
as the radial distribution of Cherenkov light or the trig- 
ger probabilities of the tele scope s can be verified experi- 
mentally (see also Hofmann 1997). Simultaneous sampling 
of the intensity of the Cherenkov light front in different 
locations emphasizes the calorimetric nature of the en- 
ergy determination, and reduces the effect of local fluc- 
tuations. Finally, given the unambiguous reconstruction 
of the shower geometry and the fact that at energies of 
one TeV virtually all events within 100 m from the cen- 
tral telescope trigger the system, and that above a few 
TeV almost all events within 200 m trigger, the effective 
detection area above 1 TeV can be basically defined by 
pure geometry, without relying on simulations. Only the 
threshold region requires a critical consideration. 

Under ideal conditions, the differential energy spec- 
trum of the incident radiation is determined as 



4>{E) 



r{E) 



V (E)A(E) 



(1) 



where r(E) is the measured rate of 7-rays of energy E 
after background subtraction, A(E) is the effective detec- 
tion area and r)(E) is the efficiency of the cuts applied 
to isolate the signal and to suppress the background. The 
effective area and the cut efficiencies are usually derived 
from Monte Carlo simulations. 

A complication arises from the finite energy resolution, 
described by the response function TZ(E' , E), the probabil- 
ity that a 7-ray of energy E is reconstructed at an energy 
E' . The measured rate is hence given by the convolution 



r(E') = I dE n{E',E) r){E)A{E)<j)(E) 



(2) 



E 

<" 




1 Replacing the 200 m by a 300 m restriction, the number 
of events reconstructed between 5 TeV and 10 TeV. and above 



Fig. 1. Effective area of the HEGRA IACT system as a 
function of energy, for vertical 7-rays. The saturation at 
10 5 m 2 at high energies reflects the cut in impact distance 
at 200 m relative to the central telescope. The uncertainty 
of the effective area in the threshold region caused by a 
5% variation of the detection threshold and by the inter- 
polation in zenith angles is shown by the hatched area. 
The corresponding systematic error on the effective area 
is - 50% at 500 GeV and - 10% at 1 TeV. 



Eq. y can no longer be trivially inverted to yield <fi(E). 
Options to find 4>(E) include the explicit deconvolution 
using a suitable algorithm, which will usually make some 
assumption concerning the smoothness of the spectrum. 
Another approach is to assume a certain functional form 
for the shape of the spectrum, and to determine free pa- 
rameters such as the flux and the spectral index from a 
fit of Eq. ^ to the data. Finally, one can absorb the effect 
of the energy smearing into a modified effective area A, 
defined such that Eq. |l| holds. The latter approach is the 
simplest, but has the disadvantage that now A depends 
on the assumed shape of the spectrum. However, with the 
< 20% energy resolution provided by the HEGRA CT sys- 
tem, shape-dependent corrections are negligible for most 
practical purposes, and results are stable after one itera- 
tion. Therefore, while both other techniques were pursued, 
the final results are based on this third method. 

The typical energy dependence of the effective area is 
shown in Fig. [|. Below 1 TeV, the effective detection area 
rises steeply with energy, and then saturates at around 
10 5 m 2 . The saturation reflects the cut on a maximum 
distance from the central telescope of 200 m. The techni- 
cal implementation of the energy reconstruction and flux 
determination is described in detail in Paper 1 . Compared 
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angle which varies during runs and from run to run. 
In order to be able to interpolate between Monte Carlo- 
generated effective areas at certain discrete zenith angles, 
a semi-empirical scaling law is exploited, which relates the 
effective areas at different zenith angles. The variation of 
zenith angles with time is accounted for through replac- 
ing Eq. [l] by the sum over all events recorded within the 
observation time T 



HE) 



1 

T 



E 



i 



r)(E,6)A(E,6) 



(•3) 



where each event is weighted with the appropriate effective 
area, given its energy and zenith angle. Note that for each 
period of a certain hardware configuration a set of effective 
areas is used which has been determined from the Monte 
Carlo simulations which model in detail the specific hard- 
ware performance, i.e. which take into account the trigger 
configuration and the mirror point spread function (see 
Paper 1). 

The key aspect in a precise and reliable determination 
of 7-ray spectra is the control of systematic errors. Sources 
of systematic errors include, e.g., 

— Systematic errors in the determination of the absolute 
energy scale. 

— Deviations from the linearity of the energy reconstruc- 
tion, caused e.g. by threshold effects at very low en- 
ergies, and possible saturation effects in the PMTs or 
the electronics at very high energies. 

— Systematic errors in the determination of the effec- 
tive area A; particularly critical is the threshold region, 
where A is a very steep function of E. 

— Systematic errors in the determination of the (energy- 
dependent) efficiency of the angular and image shape 
cuts. 

As discussed in detail in Appendix A, a non-accurate mod- 
eling of the detector response or the atmospheric trans- 
mission possibly results in 1) a shift in the energy scale of 
the reconstructed 7-ray spectrum, and 2) a distortion of 
the shape of the spectrum. Most systematic uncertainties, 
e.g. mirror reflectivities and PMT quantum efficiencies, 
exclusively contribute to an error of the energy scale. The 
calculations in Appendix A show that the curvature of 
the spectrum is reconstructed correctly provided that the 
Monte Carlo simulations accurately model the correlation 
between the detector threshold and the reconstructed en- 
ergies including the fluctuations involved. If the Monte 
Carlo description of this correlation is incorrect, the en- 
ergies reconstructed in data and the effective areas com- 
puted from Monte Carlo simulations do not match. A shift 
of the reconstructed energies by a factor (1 + e) relative 
to the effective area A used for the evaluation of the spec- 
trum results in a flux error of 



i 
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Fig. 2. Measured light yield as a function of distance to 
the shower core for 7-ray showers in the energy range of 
about 0.9 to 1.8 TeV compared with Monte Carlo simu- 
lations (hatched band). Note that the light yield assigned 
to images depends on the field of view of the camera, and 
on the definition of "image" pixels; the data should not be 
compared to "raw" simulations not including such effects. 



The flux error is potentially large in the threshold region, 
where A varies quickly with E, A cx E@ with /3 ~ 6 and 
Ac/)/<j) 6e, but is negligible at high energies, where A is 
constant, and is simply governed by the geometrical cuts 
(see Fig. 1). 

The remainder of this section is dedicated to a discus- 
sion of the various sources of systematic errors. 

4--1. Reliability of the determination of the shower energy 

A crucial input for the determination of shower energies 
is of course the expected light yield as a function of core 
distance. Since the average core distance varies signifi- 
cantly with energy, inadequacies in the assumed relation 
will not only worsen the energy resolution, but will sys- 
tematically distort the spectrum. With the redundant in- 
formation provided by a system of Chcrcnkov telescopes, 
it is possible to actually measure the light yield as a func- 
tion of the distance from the shower core and to verify the 
simulations (Aharonian et al. 1998). Briefly, the idea is to 
select showers with a fixed impact parameter relative to a 
telescope A, and with a fixed light yield in this telescope. 
This provides a sample of showers of constant energy, and 
now the light yield in other telescopes can be measured 
as a function of core distance. Fig. ^ shows the charac- 
teristic shape of the light pool for 7-ray showers, which 
is well reproduced by the simulation; this holds also for 
the variation of the sha pe wi th shower energy and zenith 
angle (Aharonian et al. |l 99S ) . 

A second key ingredient in the energy determination is 
the reconstruction of the core location. Unlike in the case 
of the angular resolution, where it is easy to show that 
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Fig. 3. (a) Difference in the x coordinate of the shower 
core as measured by independent subsystems of two tele- 
scopes for events where all four telescope triggered. Show- 
ers are selected to provide a minimum stereo angle of 
20° in each subsystem, and a maximum core distance of 
200 m from the central telescope. Points show the data, 
the dashed line the simulation, and the full line a Gaussian 
fit to the data with a width of 14 m. (b) Comparison of 
the energies measured by the two subsystems. The width 
of the Gaussian fit is 25%. Same cuts as in (a). 



the simulations, a direct check of the precision of the core 
reconstruction is not possible. However, noting that two 
telescopes suffice for a stereoscopic reconstruction, one can 
split up the four-telescope system into two systems of two 
telescopes and compare the results (Hofmann 1997). Fig- 
ure ||a illustrates the difference in core coordinates bet- 
ween the two subsystems for data and for the Monte Carlo 
simulations. As can be recognized the Monte Carlo simu- 
lations accurately predict the distribution of the distances 
between the two cores. Under the assumption that the two 
measurements are uncorrelated and that the reconstruc- 
tion accuracy achieved with two telescopes is the same for 
two telescope and four telescope events the width of the 
difference distribution should be \[2 times the resolution 
of a two-telescope system. By this means, the 2 telescope 
resolution is determined to be 14 m / \[2 w 10 m. The 
Monte Carlo simulations show, that this reconstruction 



The same technique can be used to study the energy 
resolution (Fig. |^b; see also Hofmann 1997). Also here 
the data is in excellent agreement with the simulations. 
In this case, the Monte Carlo studies predict, that the 
energies measured with the two subsystems are consider- 
ably correlated: assuming uncorrelated estimates, an en- 
ergy resolution of 25% / \/2 = 18% is inferred; the Monte 
Carlo simulations predict a true energy resolution for two 
telescope events of 23%. On the basis of simulations, fluc- 
tuations in the shower height can be identified as the ori- 
gin of this correlation. Work is ongoing to use the stereo- 
scopic determination of the shower height to improve the 
energy resolution. The excellent agreement between data 
and Monte Carlo simulations confirms that the experi- 
mental effects entering the energy determination are well 
under control, and that consequently the estimate of the 
energy resolution based on the simulations is reliable. 

4-2. The absolute energy scale 

The absolute energy calibration of IACTs is a significant 
challenge, lacking a suitable monoenergetic test beam. 
Factors entering the absolute energy calibration are the 
production of Cherenkov light in the shower, the prop- 
erties and the transparency of the atmosphere, and the 
response of the detection system involving mirror reflec- 
tivities, PMT quantum efficiencies, electronics calibration 
factors etc. So far, mainly three techniques were used to 
calibrate the HEGRA telescopes: 



The comparison between predicted and measured 
cosmic-ray detection rates. Given the integral spectral 
index of 1.7 for cosmic rays, an error e in the energy 
scale results in an error of 1.7e in the rate above a given 
threshold. Apart from the slightly different longitudi- 
nal evolution of hadronic and of 7-ray induced show- 
ers, this test checks all factors entering the calibration. 
The Monte Carlo simulations reproduce the measured 
cosm ic-ray trigger rates within 10% (Konopelko et al. 
1999 ). The I ACT system was furthermore used to de- 



termine the flux of cosmic-ray protons i n the 1.3 to 
10 TeV energy range (Aharonian et al. 1999b| ). The 
measured flux of protons 



dN p /dE = (0.11 ± 0.02 stat ± 0.05 sys )x 



(5) 



E 



-2.72±0.02 stat ±0.15 sv 



'TeV" 



excellently agrees with a fit to the combined data of 
all other experiments, indicating that systematics are 
well under control and that the assigned systematic er- 
rors - partially related to the energy scale - are rather 
conservative. 

Given the measured characteristics of the telescope 
components such as mirrors or PMTs, the sensitivity 
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Fig. 4. Mean reconstructed shower energy as a function 
of the true energy for showers incident under 20° zenith 
angle, also showing the rms errors of the energy recon- 
struction. 



Using a distant, calibrated, pulsed light source, an 
overall calibration of the response of the telescope and 
its readout electronics co uld b e achieved, with a pre- 
cision of 10% (Frafi et al. |1997|) . 



The last three techniques have to rely on the Monte 
Carlo simulations of the shower and of the atmospheric 
transparency. While Monte Carlo simulations have con- 
verged and different codes produce consistent results, de- 
tails of the atmospheric model and assumptions concern- 
ing aerosol densities can change the Cherenkov light yield 
on the ground by about 8% (Hemberger 1998). 

Within their errors, all calibration techniques are con- 
sistent. Overall, we believe that a 15% systematic error 
on the energy scale is conservative, given the state of sim- 
ulations and understanding of the instrument. Based on 
the comparison with cosmic-ray rates, one would conclude 
that the actual calibration uncertainty is below 10%. 

4-3. The threshold region and associated uncertainties 

As discussed earlier in detail, the threshold region is 
very susceptible to systematic errors. In the sub-threshold 
regime events trigger only because of upward fluctuations 
in the light yield, and energy estimates tend to be biased 
towards larger values. Fig. ^ shows the mean reconstructed 
energy as a function of the true energy for a sample of sim- 
ulated events at typical zenith angles. Note that for the 
1997 data set we have a noticeable number of detected 7- 
rays events below 500 GeV. However in this energy region 
the bias is so strong that a reliable correction is no longer 
possible. Therefore spectra will only be quoted above this 
energy. 




Fig. 5. Size-distributions of 7-ray induced air showers for 
data and simulations for two of the four telescopes. The 
data is from a period of a certain trigger configuration, 
namely "Data-period F of Paper 1, and the Monte Carlo 
simulations used the telescope specific single pixel trigger 
probability as function of signal amplitude derived from 
this data (Mkn 501 data, background subtracted, Monte 
Carlo weighted according to the results in Sect. 5). 



the simulation great emphasis was placed on the correct 
description of the pixel trigger probabilities as a function 
of signal amplitude. For each trigger configuration and for 
each telescope this dependence has been derived from air 
shower data using the recorded information about which 
pixels of a triggered telescope surpassed the discriminator 
threshold and which pixels did not. As discussed above 
and in Appendix A, the onset of the size-distribution pro- 
vides a sensitive test of the quality of the simulation. Fig. |^ 
shows the measured and simulated size-distribution for 
two telescopes for a certain trigger configuration ("Data- 
period I" of Paper 1). For our current best simulations, fits 
to the rising edges of the size-distributions indicate for the 
individual telescopes and the different trigger configura- 
tions deviations of the size scale between data and Monte 
Carlo on the 5% level. For conservatively estimating the 
systematic error on the shape of the Mkn 501 spectrum, 
we allow for a ±5% correlated shift of the thresholds of all 
four telescopes. The resulting uncertainties in the effective 
detection area A are indicated in Fig. EL as expected, A 



8 



F. Aharonian et al.: The 1997 Time Averaged TeV Spectrum of Mkn 501 



Random threshold variations (Gaussian-distributed 
with a width of 15%) between individual pixels were found 
to be of relatively small influence compared to the system- 
atic threshold shifts. 

As another source of systematic errors of special im- 
portance in the threshold region, the accuracy of the en- 
ergy reconstruction for zenith angles between the discrete 
simulated zenith angles {9 = 0°, 20°, 30°, 45°) has been 
considered. Imperfections in the scaling law used to re- 
late the Cherenkov light yield at different zenith angles 9 
could result in a systematic shift of the reconstructed en- 
ergy at intermediate values of 9. To test the description, 
the energy of the Monte Carlo showers with 9 = 30° was 
determined with the light yield tables of the 0°- and 45°- 
showers and the result was compared to the result based 
on using the 30° light yield table. The systematic shift 
of the reconstructed energies was about 5% below 1 TeV 
and 2% above 1 TeV; based on this and other studies we 
believe that using the full set of simulations, systematic 
shifts in the reconstructed energy are below 5% and 2%, 
below and above 1 TeV respectively. 

The modified effective area used for the determination 
of the spectrum slightly depends on the assumed source 
spectrum. At the lowest and at the highest energies the 
source spectrum can not be determined with high statisti- 
cal accuracy. For energies below 1 TeV we conservatively 
estimate the corresponding uncertainty in the modified ef- 
fective area, by varying the spectral index of an assumed 
source spectrum dN/dEcxE~ Q from a = 1.5 to a = 3. At 
the highest energies above 15 TeV we vary the assumed 
source spectrum from a broken power law to a power law 
with an exponential cutoff, both specified by fits to the 
data. In addition we explore the statistical significance of 
the 7-ray excess at the highest energies by a dedicated 
X 2 -analysis (see Sect. |J). 

In the threshold region of the detector the total sys- 
tematic error is dominated by the systematic shift of the 
telescope thresholds and by the possible systematic shift in 
the reconstructed energy due to the zenith angle interpola- 
tion. The total systematic error is computed by summing 
up the individual relative contributions in quadrature. 

4-4- Saturation effects at high energies 

The PMTs and readout chains of the HEGRA system tele- 
scopes provide a linear response up to amplitudes of about 
200 photoelectrons. At higher intensities, nonlinearities of 
the PMT become noticeable, and also the 8-bit Flash- 
ADC saturates. With typical signals in the highest pixels 
of 25 photoelectrons per TeV 7-ray energy, the effects be- 
come important for energies around 10 TeV and above. 

Saturation of the Flash- ADC can be partly recovered 
by using the recorded length of the pulse to estimate 
its amplitude, effectively providing a logarithmic charac- 
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Fig. 6. The mean signal in the peak pixel, hoti, as a func- 
tion of the image size, for data (full points) and simula- 
tions (open points). 



and a correction is applied. These nonlinearities and the 
Flash- ADC saturation are included in the simulations. As 
a very sensitive quantity to test the handling of satura- 
tion characteristics, the dependence of the pulse height in 
the peak pixel on the image size emerged (Fig. |^). Data 
and simulations are in very good agreement up to pixel 
amplitudes exceeding 10 3 photoelectrons, equivalent to 
>50 TeV 7-ray showers. Older versions of the simulations, 
which did not properly account for PMT/preamp nonlin- 
earities, showed marked deviations for size-values above 
10 3 . 

Independent tests of saturation and saturation correc- 
tions were provided by omitting, both in the data and 
in the simulation, the highest pixels in the image, and 
by comparing event samples in different ranges of core 
distance and zenith angle. The saturation effect is non- 
negligible only at very high energies, namely E > 15 TeV. 
However, all tests show that given the quality of our cur- 
rent simulations, systematic errors induced by saturation 
effects even in this energy region are yet small compared 
to the statistical errors. 

4-5. Efficiency of cuts 

Among the sources of systematic errors, the influence of 
cuts is less critical. The large flux of 7-rays from Mrk 501 
combined with the excellent background rejection of the 
IACT system allows to detect the signal essentially with- 
out cuts. In the analysis, one can afford to apply only 
rather loose cuts, which keep over 80% of all 7-rays; only 
at the lowest energies, a slightly larger fraction of events is 
rejected. Since only a small fraction of the signal is cut, the 
uncertainty in the cut efficiency is a priori small; in addi- 
tion, the efficiency can be verified experimentally by com- 
paring with the signal before cuts, see, e.g., Fig. IH The cut 
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Fig. 7. Efficiency of the cut on the shower direction rela- 
tive to the Mrk 501 location after the software threshold 
of at least 40 photoelectrons in two or more telescopes (a) , 
and of the shape cuts used to enhance 7-rays (b), as a func- 
tion of energy. Full points show the measured efficiencies, 
open points the results of the Monte Carlo simulations. 
Efficiencies determined from data can be larger than one 
due to background fluctuations. 



simulations. We conservatively estimate the systematic er- 
ror the cut efficiencies, to be 10% at 500 GeV decreasing 
to a constant value of 5% at 2 TeV and rising above 10 
TeV to 15% at 30 TeV. At low energies the acceptances 
are corrected according to the measured efficiencies. 

4-6. Other systematic errors and tests 

The precision in the determination of the effective areas 
is partly limited by Monte Carlo statistics. The universal 
scaling law used to relate Monte Carlo generated effec- 
tive areas at different zenith angles involves a rescaling of 
shower energies. Statistical fluctuations in a Monte Carlo 
sample at a given energy and angle will hence influence 
the area over a range of energies. Because of the resulting 
slight correlation, Monte Carlo statistics is in the following 
included in the systematic errors. 

To test for systematic errors, the data sample was 
split up into subsamples with complementary systematic 
effects, and spectra obtained for these subsamples were 
compared. Typical subsamples include 

— Events where 2, 3 or 4 telescopes are used in the re- 
construction. 




10 

Energy [ TeV ] 

Fig. 8. Ratio of the spectra computed with the events 
with small (0 to 120 m) and with large (120 m to 200 
m) impact distances relative to the central telescope. A 
fit to a constant gives a ratio of 1.04 ± 0.02 with a x 2 of 
23.1 for 12 degrees of freedom. The fit is restricted to the 
energy region of reasonable systematic errors, i.e. with an 
effective detection area for both event samples larger than 
~ 10 4 m 2 . 

or a signal in the two peak pixels of at least 30 photo- 
electrons. 

— Events with showers in a certain distance range from 
the center of the system (CT3), e.g. 0-120 m compared 
to 120-200 m. This comparison tests systematics in 
the light-distance relation as well as the correction of 
nonlinearities in the telescope response. 

— Events where all pixels are below the threshold for non- 
linearities. 

— Events in different zenith angle ranges. The compar- 
ison of these spectra provides a sensitive test of the 
entire machinery, and also of nonlinearities, where the 
data at larger angles should be less susceptible because 
of the smaller size at a given energy. 

For each of the subsamples, the effective area and cut effi- 
ciencies were determined, and a flux was calculated. In all 
cases, deviations between subsample spectra were insignif- 
icant, or well within the range of systematic errors. Among 
the variables studied, the most significant indication of 
remaining systematic effects is seen in the comparison of 
different ranges in shower impact parameter relative to 
the central telescope. The ratio of the spectra determined 
with the data of small (<120 m) and large (between 120 
m and 200 m) impact distances is shown in Fig. ||. A fit 
to a constant gives a mean ratio of 1.04 ± 0.02 with a % 2 - 
value of 23.1 for 12 degrees of freedom, corresponding to 
a chance probability for larger deviations of 5%. 

5. Experimental results 

In this section we first present the 1997 Mkn 501 time 
averaged energy spectrum. As discussed already in the in- 
troduction the derivation of a time averaged spectrum is 
meaningful since the changes in the spectral shape dur- 
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Fig. 9. Time-averaged energy spectrum of Mrk 501 for 
the 1997 observation period. Vertical errors bars indicate 
statistical errors. The hatched area gives the estimated 
systematic errors, except the 15% uncertainty on the ab- 
solute energy scale. The lines shows the fit discussed in 
the text. 



cally between 0.1 and 0.3 in the diurnal spectral indices. 
Moreover, as described in Paper 1, dividing the data into 
groups according to the absolute flux level or according to 
the rising or falling behavior of the source activity yielded 
mean spectra which did not differ significantly from each 
other in the one to ten TeV energy range. The weakness 
of the correlation between the absolute flux and the spec- 
tral shape will further be substantiated below over the 
energy region from 500 GeV to 15 TeV. Nevertheless, the 
importance of the spectral constancy should not be over- 
estimated. If the spectral variability is not tightly corre- 
lated with the absolute flux, diurnal spectral variability 
characterized by a change of the spectral index at sev- 
eral TeV by approximately ±0.1 is surely consistent with 
the HEGRA data. The time-averaged energy spectrum is 
shown in Fig. |9|. For the determination of the spectrum 
also at energies below 800 GeV, only the data from zenith 
angles smaller 30° have been used (80 h observation time). 
The measurements extend from 500 GeV to 24 TeV. The 
hatched region in Fig. || ff. gives our estimate of the sys- 
tematic errors on the shape of the spectrum, except the 
15% uncertainty on the absolute energy scale. The spec- 
trum shows a gradual steepening over the entire energy 
range. A fit of the data from 500 GeV to 24 TeV with a 
power law model with an exponential cut off gives: 



Table 1. The time-averaged differential spectrum of Mkn 
501 
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a energy in TeV 
-l 



b in (cm" 2 s ic v 
statistical error in (cm -2 s" 1 TeV" 1 ) 
systematic error on the shape of the spectrum in 
(cm" 



s -i TeV" ) 
upper limits in (cm" 



TeV" 



at 2 a confidence level 



a - 



(10.8 ±0.2 stat ±2.1. 



1.92 ± 0.03 s 



± 0.20 s 



10" 

, and Eq 



^TeV" 
(6.2 



± 



0.4 sta t (—1.5 + 2.9) sys ) TeV. The systematic errors on 
the fit parameters result from worst case assumptions con- 
cerning the systematic errors of the data points, and their 
correlations and include the error caused by the 15% un- 
certainty in the energy scale. The errors on the fit param- 
eters, especially on a and Eq, are strongly correlated. The 
variation of only one of the parameters within the quoted 
error range yields spectra which are inconsistent with the 
measured spectrum. The data points and their errors are 
summarized in Table 1. 

In the highest energy bin (19 TeV to 24 TeV) 40 excess 
events are found above a background of 13 events, corre- 
sponding to a nominal significance of S = (N on - N g) / 
VNon + N ff of 3.7 a. However, due to the steep spectrum 
in this energy range, a part of these events may represent 
a spill-over from lower energies. To provide an absolutely 
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Fig. 10. The spectral energy distribution E 2 dN/dE, for 
the data set of low zenith angles (9 < 30° , full circles) and 
for the data set of large zenith angles (8 between 30° and 
45°, 32 h observation time, open symbols). Since the obser- 
vation periods do not overlap for the variable source, the 
spectra are normalized at the energy 2 TeV. The hatched 
band indicates the systematic error on the shape of the 
spectrum for the low zenith angle data. The systematic 
error on the high zenith angle spectrum at energy E ap- 
proximately equals the systematic error on the low zenith 
angle spectrum at energy E/2. 



a sharp cutoff at E = E cut : dN/dE = N a (E/l TeV)" Q 
exp {—E/Eq) Q(E cut — E). The best fit is achieved with 
E cut = 28 TeV; the 2a lower limit is E cut = 16 TeV. 

Fig. |l0| illustrates the spectral energy distribution, 
E 2 dN/dE(E) as determined from the small zenith an- 
gle data (<30°, energy threshold 500 GeV) and the large 
zenith angle data (30° to 45°, 32 h observation time, en- 
ergy threshold 1 TeV). Note that the large zenith angle 
data has mainly been acquired during the second half of 
the 1997 data taking period. Nevertheless the shape of 
both spectra agrees within the statistical and systematic 
errors. The combined small and large zenith angle data set 
yields the same lower limit on E cut of 16 TeV as derived 
from the small zenith angle data alone. It can be recog- 
nized that the spectral energy distribution is essentially 
flat from 500 GeV up to ~ 2 TeV. 

Figure [ll] (upper panel) shows the spectral energy 
distribution for the overall data sample and for peri- 
ods of high and low flux separately: dN/dE (2 TeV) 
determined on diurnal basis above 3 and below 1.6 
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Fig. 11. The upper panel illustrates the spectral energy 
distribution E 2 dN/dE, for the full data set (full circles), 
for periods of low flux (open circles), and for periods of 
high flux (triangles) ( dN/dE(2 TeV) above 30 and below 
16 times 10 _12 cm _1 s _1 TeV _1 ). Only the statistical er- 
rors are given here; the systematic errors enter the three 
spectra in the same way and can be neglected compar- 
ing the three spectra. The dashed lines indicate the shape 
of the mean spectrum (fit from Eq. ^|) overlaid over all 
three spectra to simplify the comparison of the shape of 
the three spectra. In the lower panel the ratio of the low 
flux spectrum divided by the high flux spectrum is shown. 
The dashed line gives the fit to a constant. The x 2 -value 
is 12.3 for 15 degrees of freedom. 



errors, as shown by the ratio of both spectra, presented 
in Fig. [ll] (lower panel). The systematic error is to good 
approximation the same for both data samples and can- 
cels out in the ratio. The result thus confirms our previous 
conclusion about the flux-independence of the spectrum of 
Mkn 501 in 1997 between 1 and 10 TeV (Paper 1). Now 
the statement is extended to the broader energy region, 
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Fig. 12. The Time-averaged spectrum of Mrk 501 dur- 
ing 1997, compared with published results from other ex- 
periments (Aharonian et al. 1999c, Hayashida et al. 1998 
Samuelson et al. 1998, Barrau et al. 1998). Since the obser- 



vation periods do not completely overlap for the variable 
source, the spectra are normalized at the energy 2 TeV. 
For the HEGRA system the hatched area shows the sys- 
tematic errors on the shape of the spectrum as described 
in the text. For the other experiments only statistical er- 
rors are shown. Only data points with a signal to noise 
ratio larger than one have been used. 



accuracy, e.g. a power law fit in the energy region from 
1 TeV to 5 TeV gives a differential index of -2.23 ±0.04 stat 
and -2.26 ±0.06 s t a t for the high and the low flux spectrum 
respectively. In the narrow energy range from 500 GeV to 
1 TeV the statistical uncertainty on the spectral index 
is considerably larger, we compute 0.2 for the high flux 
sample and 0.4 for the low flux sample. Therefore, our 
1997 Mkn 501 data would not contradict a correlation of 
emission strength and spectral shape below 1 TeV as ten- 
tatively reported by the CAT-group (Djannati-Atai et al. 
19991) . 



For completeness, the HEGRA IACT system data are 
plotted in Fi g. P j ointly with the HEGRA CTri (Aha- 



ronian et al. |1999cj ), the CAT (B arrau et al. |1998| ), the 



Telescope Array (Hayashida et al 
pie (Samuelson et 



al. 



1998), and the Whip- 



1998) results concerning the Mrk 



501 energy spectrum during the 1997 outburst. Generally 
a good agreement can be recognized in the overlapping 
energy regions, except for a steeper Telescope Array spec- 
trum. 



6. Discussion 

The observations of Mkn 501 by the HEGRA IACT sys- 
tem during the long outburst in 1997 convincingly demon- 
strate for the first time that the energy spectrum of the 
source extends well beyond 10 TeV. We believe that this 
very fact, together with the discovery of a time- and flux- 
independent stable spectral shape of the TeV radiation 
will have considerable impact on our understanding of the 
nonthermal processes in AGN jets. 

This discussion is not an attempt at detailed modeling 
of the result presented in the previous section; this will 
be done elsewhere. Neither shall we systematically invoke 
multi-wavelength observations. Our purpose is rather to 
point out the multiple facets of the 7-ray phenomenon on 
its own. They stem from the fact that the emission prob- 
ably originates from a population of accelerated particles 
within a spatially confined relativistic jet, specifically ori- 
ented towards the observer, and that subsequently the ra- 
diation must propagate through the diffuse extragalactic 
background radiation field (DEBRA) before it reaches us. 
Each of these circumstances can influence the character- 
istics of the emitted spectrum, and we shall address them 
in turn below. 

6.1. Production and absorption 0/ TeV photons in the jet 

The enormous apparent VHE 7-ray luminosity of the 
source, reaching ~ 10 45 erg/s during the strongest flares 
which typically last At = 1 day or less, implies that 
the 7-rays are most probably produced in a relativistic, 
small-scale (sub-parsec) jet which is directed along the 
observer's line of sight. The determining quantity is the 
so-called Doppler factor 8- 3 = [(1 + z)T- } {\ — /3j cos 6)] -1 , 
where z is the redshift of a source moving with velocity 
/3jc and Lorentz factor T s = (1 — /3?) -1 / 2 along the jet axis 
that makes an angle O with the direction to the observer. 
Moreover, the assumption of relativistic bulk motion ap- 
pears to be unavoidable in order to overcome the problem 
of severe 7 — 7 absorption by pair production on low- 
frequency photons inside the source (see e.g. Dermer & 
Schlickciscr 1994). Indeed, assuming that the 7-radiation 



is emitted isotropically in the frame of a relativistically 
moving source, the optical depth at the observed 7-ray 
energy is easily estimated as 



'77 



frd 2 a T 5r b E 
8m2 c 6 At 



0.065/_xoAtTi5- 6 if- 2 ^ 



(7) 



Here S w = <5j/10, At day = At/1 day, E TcV = E/lTeV, 
with E the energy of the gamma-ray in the laboratory 
frame, and /_io = /r/10~ 10 erg/cm 2 s is the observed en- 
ergy flux at hv ~ 100(5 2 EZ e y eV which for an order of 
magnitude estimate is assumed to be constant at the op- 
tical to UV wavelengths that predominantly contribute to 
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normalized to the value of the Hubble constant Hqo = 
i?o/60km/sMpc. Assuming now that the observed opti- 
cal/U V flux of Mkn 501, f- W ~ 0.5 (see e.g. Pian et al. 
1998| ), is produced in the jet, the absorption of 20 TeV 
7-rays becomes negligible only when <5j > 10; already 
for <5j = 8 internal absorption would be catastrophic for 
A^day = 1, with e _T ^ < 10~ 2 . One may interpret the 
detected VHE spectrum of Mkn 501 given by Eq. || as a 
power-law production spectrum, modified by internal 7 — 7 
extinction with optical depth r = i?Tcv/6.2. This would 
give an accurate determination of the jet's Doppler factor 
taking into account the very weak dependence of Sj on all 

relevant parameters, namely Sj = 8.5 /lio^da/ 6 ^-60^- 
Since there could be a number of other reasons for the 
steepening of the TeV spectrum, that estimate can only 
be considered as a robust lower limit on <5j . 

The strong steepening of the observed spectrum of 
Mkn 501 above several TeV could also be attributed, for 
example, to an exponential cutoff in the spectrum of ac- 
celerated particles. These could either be protons produc- 
ing 7-rays through inelastic p — p interactions and sub- 
sequent 7r°-decay , or electrons producing 7-rays via in- 
verse Compton (IC) scattering. In IC models an additional 
steepening of the 7-ray spectrum is naturally expected 
due to the production of TeV 7-rays in the Klein-Nishina 
regime if the energy losses of electrons are dominated by 
synchrotron radiation in the jet's magnetic field. And fi- 
nally, the exponential cutoff in the observed TeV spec- 
trum could be caused by 7-7 absorption of TeV photons 
in a possible dust torus surrounding the AGN (Protheroe 
& Biermann 1996|) and in the ex tragalactic diffuse back- 
ground radiation (Nikishov 1962 Gould & Schreder 1965| , 



1996 



Jelly 1965, Stecker et al. 1992) 



Our current poor knowledge about the distortion of the 
source spectrum caused by internal and intergalactic ab- 
sorption does not allow us to distinguish between hadronic 
and leptonic source models on the basis of their predic- 
tions concerning the TeV energy spectra. Fortunately the 
temporal characteristics are to a large extent free from 
these uncertainties. Thus we believe that real progress in 
this area can only be achieved by the analysis of both the 
spectral and temporal characteristics of X-ray and TeV 7- 
ray emissions obtained during multiwavelength campaigns 
that investigate several X-ray selected BL Lac objects in 
different states of activity, and located at different dis- 
tances within several 100 Mpc. Notwithstanding this be- 
lief, we show here that the high-quality HEGRA spectrum 
of Mkn 501 alone allows us to make quite a few interest- 
ing inferences about the 7-ray production and absorption 
mechanisms. 



6.1.1. IC models of the gamma ray emission. 



1997 



length observations of Mkn 421 (Takahashi et al. 
Buckley et al. |l996j ) and Mkn 501 (Catanese et al 
Pian et al. [1998 , Paper 1), support the hypothesis of both 
emission components originating in relativistic jets due to 
synchrotron/IC radiation of the same pop ulation of di- 
rectly accele rated electrons (Ghisellini et al 
1996| ; Inoue & Takahara 



Mars cher 
Kirk Il997 



1996 



1996: Bloom 



19971) 



Mastichiadis 
One of the 



Bcdnarek & Protheroe 
distinctive features of leptonic models is that they allow 
significant temporal and spectral variations of TeV radi- 
ation. The stable shape of the TeV spectrum of Mkn 501 
during the 1997 outburst does not contradict these mod- 
els. It rather requires them to have two important features. 

First of all the form of the spectrum of accelerated 
electrons should be essentially stable in time and be inde- 
pendent of the strength of the flare up to electron energies 
responsible for the production of the highest observed 7- 
ray energies E > 16TeV. For a typical Doppler factor 
<5j ~ 10 this implies relatively modest electron energies 
E c < 10 TeV in the jet frame, assuming that the Comp- 
ton scattering at the highest energies takes place in the 
Klein-Nishina limit. In this limit a significant fraction of 
the electron energy goes to the upscattered photon, i. e. 
/•• - »■ I - ■ 

For low values of the magnetic field in the emitting 
plasma, i.e. B < 0.01 G, the X-ray spectrum could be 
more sensitive to accelerated electrons with energies above 
10 TeV. The typical observed energy of X-rays produced 
by electrons of energy E e in the jet frame is 



E x ^ 20(B/0.1 G) (E c /1 TeV) 2 6 10 keV. 



(8) 



The BeppoSAX observations of Mkn 501 in April 1997 
showed that the X-ray spectrum becomes very hard dur- 
ing strong flares. This is interpreted as a shift of the syn- 
chrotron peak to energies in excess of 100 keV (Pian et 
al. 1995). Formally this effect could be explained by a sig- 
nificant increase in each of the three parameters which 
determine the position of the synchrotron peak, i.e. the 
maximum electron energy E c _ max , the magnetic field B, 
and the jet Doppler factor 8 y 

The rather stable energy spectrum of TeV radiation 
implies that the spectrum of the parent electrons does 
not significantly vary during the HEGRA observations, 
the latter performed with typical integration times be- 
tween one and two hours. The condition of a constant 
acceleration spectrum does not yet guarantee a stable 
energy spectrum of TeV radiation. Therefore we need a 
second condition, namely to assume very effective radia- 
tive (synchrotron and IC) cooling of electrons, sufficiently 
fast to establish an equilibrium electron spectrum within 
At* = 10 <5io h. The radiative cooling time is 



^rad = (o CT T C W E c /m Q C 2 ) 1 



(9) 
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where wq = B 2 /8ir + w T is the total energy density of 
magnetic and photon fields. Thus, for a jet magnetic field 
of about 0.1 G and a comparable low- frequency photon 
density (w 4- 10 -4 erg/cm 3 ) a radiative cooling time of less 
than 5 hours (in the jet frame) could be easily achieved. 

6.1.2. 7T° origin of gamma rays 

The lack of correlation between spectral shape and ab- 
solute flux, as well as the very fact that 7-rays with en- 
ergy > 16 TeV are observed, could also be explained, per- 
haps even in a more natural way, by the assumption of 
a '7r°-decay' origin of the 7 -radiation. Yet the efficiency 
of this mechanism in the jet appears to be too low to ex- 
plain the observed time variability and the high fluxes of 
the TeV radiation. This is due to the low density 7ih of 
the thermal electron-proton plasma in the jet. The prob- 
lem of variability could be at least in principle overcome 
by invoking adiabatic losses caused by relativistic expan- 
sion of the emitting "blob". However, this assumption 
implies very inefficient 7-ray production with a luminos- 



T t If* 



" 10 erg/cm 2 s is the total observed radia- 
Pian et 



al. 

-10 



1998). Then the observed 



lty L 7 

tivistic protons, i a d > R/c ~ At<5j is the adiabatic cool- 
ing time, and t pp ~ 5 • 10 15 (nii/lcm^ 3 )^ 1 s denotes the 
characteristic emission time scale of 7r°-decay 7-rays. We 
shall assume here that the proton luminosity L p should 
not exceed the total power of the central engine, roughly 
the Eddington luminosity L E = 1.3 • 10 45 (M/Mq) erg/s 
of a supermassive Black Hole of mass M — 10 7 M Q , or 
more empirically, the apparent (47r) total luminosity of 
the source which is L to t = 4 • ttcP ■ /tot ~ 10 45 erg/s, where 
/tot ~ few x 10"" 
tive flux (see e.g. 

TeV-flux of about 2.5 • 10~ 10 erg/s requires a lower limit 
n a > 10 6 cm~ 3 on the density of the thermal plasma in 
the jet. This makes the relativistically moving 'blob' very 
heavy (M ~ 0.05M Q ) with an unacceptably large kinetic 
energy E kin = Mc 2 T } ~ 10 54 Ai 3 ay <5 4 (n/10 6 cm~ 3 ) erg. 

We would like to emphasize that these arguments hold 
against the 7r°-origin of 7-rays produced in a small-scale 
jet; they do not in general exclude hadronic models. In 
particular, scenarios like the one assuming 7-radiation 
produced by gas clouds that move across the jet (Bednarek 
& Protheroe (1997p , Dar & Laor [1997D remain an attrac- 
tive possibility for hadronic models. They d o not exclude 
either a "proton blazar" model (Mannheim 1993 ). It im- 
plies a secondary origin of the relativistic electrons that 
are the result of electromagnetic cascade, triggered by 
photo-meson processes involving extremel y hig h energy 
protons in a hadronic jet (sec Mannheim 1998 and ref- 
erences therein). 

6.2. Intergalactic extinction 



astrophysically significant (see e.g. Stanev & Franceschini 
19981, Funk et al. 11998;, Biller et al. 11998;, Stecker & de 



Jager |1998| , Biller |1998| , Primack et al. |1998| , Stecker p.998| ) 



after the discovery of TeV radiation from Mkn 421 and 
Mkn 501 up to energies of 10 TeV, as reported by the 
Whipple (Zweerink et al. |l997|) , and HEGRA (Aharonian 
et al. 1997a), CAT (Djannati-Atai et al. |l999| ), and Tele- 



scope Array (Hayashida et al. 1998 ) groups 



6.2.1. Gamma ray absorption 

If we ignore the appearance of second generation 7-rays 
(see Sect. 6.2.2), then extinction is reduced to a simple 
absorption effect which can be described by a single ab- 
sorption optical depth r. 

The optical depth r of the intergalactic medium for a 
7-ray photon of energy E, emitted from a source at the 
distance d = cz/Hq, can be expressed for small redshifts 
z<C 1 in a convenient approximate form using a quantity 
t = t'£, where 



where L p is the luminosity in rela- T '{E) 



0.08 



a rc(e m ) 



(10) 



10" 3 eV/cm 3 ) 



)( 



0.034 



E Te vH. 



1,11 



iV60km/sMpc, E Tc v = E/lTeV, e m = 
li^TeV eV, and £ being a correction factor 



with Hqq 

Aml^/E 

which accounts for the specific form of the differential DE- 
BRA photon number density n(e); the background pho- 
ton energy is denoted by e. This expression is based on 
the narrowness of the 77 — > e + e~ cross-section cr 77 as a 
function of (e/e m ) which peaks at e/e m ~ 1 in an isotropic 
field of background photons (Hcrtcrich 1974). Thus for 
a large class of broad DEBRA spectra n(e) the optical 
depth is essentially caused by background photons with 
energy centered around e m . For (broad) power-law spec- 
tra, n(e) = hq6 1 , the optical depth can be calculated 
analytically as t(E) — 77(7) • 4 7 (crT/4)e m n(e m ) d, where 
77(7) = 7/67" 5 / 3 (l + 7)- 1 (Svensson |1987| ). Thus for rel- 
atively flat power-law spectra with 1 < 7 < 2.5 we ob- 
tain £ = 4 7 ?7(7) ~ 2, i.e. approximately half of r 77 is 
contributed by background photons with e in the interval 
given by e m ± 1/2 

We note that t(E) cx E 1 ^ 1 for a power-law spectrum 
of DEBRA . For example, within the 'valley' of the en- 
ergy density at mid infrared wavelengths from several /im 
to several tens of /im, where the energy density is expected 
to be more or less constant (i.e. e 2 n(e) = const), the in- 
tergalactic extinction of 7-rays is largest at the highest 
observed energies. In particular, according to Eq. 



10 



even 



at a very low and probably unrealistic level of the DE- 
BRA intensity of e 2 n(e) = 10~ 4 eV/cm 3 at A — 30 /im, 
approximately 35 percent of the 25 TeV 7-rays emitted by 
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Fig. 13. The energy fluxes of DEBRA for pure power-law 
differential spectra with 7 = 1 (curve 1), 7 = 2 (curve 2) 
and 7 = 3 (curve 3). Curve 4 is the sum of spectra 1 and 
3. The absolute flux normalizations have been determined 
from the condition of t 77 = 3 (for curve 1), and from the 
condition of the maximum possible flux of DEBRA which 
still "reproduce" reasonable 7-ray source spectra shown 
in Fig. [l4| (for curves 2 and 3). The horizontal bars cor- 
respond to the upper limits on DEBRA fluxes obtained 
using a method similar to the one suggested by Biller et 
al. (1998). The curve marked as "MBR" correspond to the 
density of the 2.7 K MBR. The tentative flux measurement 
at 3.5 /zm is taken from Dwek & Arendt (1998), and the 
flux estimates based on the ISO survey at 6 and 15 /im 
are from Stanev & Franceschini (1998). The other mea- 
sured fluxes and the upper/lower limit estimates of the 
DEBRA are taken from the recent compilation by Dwek 
et al. (1998), i.e., the upper limit at 2.2 /xm is from Hauser 
et al. (1998), the lower limit at 2.2 fim is from Gardner 
et al. (1997), a nd th e UV to optical detections are from 
Pozzetti et al. fll998|). 



contains important information about the DEBRA, at 
least at wavelengths A > 10/im. Moreover, it is quite pos- 
sible that a non-negligible intergalactic extinction takes 
place also at low 7-ray energies due to interactions with 
near infrared (NIR) background photons. Interpreting the 
power-law shape of the spectrum at 2 TeV as an indication 
for weak extinction by, say, a factor less than or equal to 2, 
and ignoring the contributions from all other wavelengths 
beyond the interval 0.5 ± 0.25 eV (A ~ 2.5lQg 5 /xm) from 
Eq. [l0[ one obtains e 2 n(e) < 4 x 10~ 3 H 60 eV/cm 3 , not far 
from the flux of DEBRA experimentally inferred (Dwek 
et al. 1998 , de Jager & Dwek 1998 ) and theoret ically ex- 
pected (Malkan & Stecker |l998j Primack et al. |l998|) at 
these wavelengths. 
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Fig. 14. The source spectra of Mkn 501 reconstructed for 
different models of DEBRA. The heavy dots correspond 
to the measured spectrum of HEGRA approximated by 
Eq. 6 with 7 = 1.92, and E = 6.2 TeV. The curves 1, 2, 
and 3 correspond to the power-law DEBRA spectra shown 
in Fig. |l3| by curves 1,2 and 3, respectively. The vertical 
line at 16 TeV indicates the edge of the 7-ray spectrum 
of Mkn 501 measured by HEGRA. 



model-independent upper limit (e.g. Weekes et al. 1997). 
Obviously for any quantitative estimate of the DEBRA 
one needs to know the intrinsic 7-ray spectrum, and this 
can not be obtained from 7-ray observations alone. As al- 
ready emphasized above, a multi-wavelength approach is 
indispensable, for example in the form of detailed model- 
ing of the entire or at least a large wave-length range of the 
nonthermal spectrum. To be specific, one avenue would be 
modeling the nonthermal X-ray and 7-ray spectra in the 
framework of synchrotron-inverse Compton models, based 
on simultaneous multi-wavelength observations of X-ray 
selected BL Lac objects (Coppi & Aharonian 1999). 

An illustration for the need of more than 7-ray data 
alone is given by a DEBRA spectrum n(e) oc e _1 . It does 
not change the spectral shape of 7-rays at all (grey opac- 
ity), although formally the extinction could be arbitrarily 
large. The curves marked as "1" in Fig. [l3] and Fig. [l4| 
may serve as an example. In the case of Mkn 501 this am- 
biguity can be significantly reduced by rather general ar- 
guments regarding 7-ray energetics. Requiring again that 
the 7-ray luminosity should not exceed the total apparent 
luminosity of the source, L to t ~ 10 45 erg/s, we must have 
e T < 10 4 ((5j/10) 4 . For 7-ray production in the jet with 
<5j ~ 10, this implies r < 10. In fact, already a value of 
r ~ 3 (corresponding to curve 1 of Fig. [ll]) creates uncom- 
fortable conditions for the majority of realistic models of 
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chrotron and inverse Compton processes. Indeed, r 77 ~ 3 
implies that the 7-ray luminosity of the source correspond- 
ing to the "reconstructed spectrum" (curve 1 in Fig. |l4| ) 
exceeds the luminosity of the source in all other wave- 
lengths by an order of magnitude which hardly could be 
accepted for any realistic combination of parameters char- 
acterizing the jet. 

Accepting the current lack of reliable knowledge of 
the 7-ray source spectrum, it is nevertheless worthwhile 
to derive upper limits on DEBRA by formulating differ- 
ent a priori, but astrophysically meaningful requirements 
on the spectrum and the 7-ray luminosity of the source. 
A possible criterion, for example, could be that within 
any reasonable source model the intrinsic spectrum of 7- 
rays, Jq(E) = J Q b s exp(r), should not contain a feature 
which exponentially rises with energy at any observed 7- 
ray energy. In practice this implies that the 'reconstruc- 
tion' factor exp(r) should not significantly exceed the 
exponential term of the observed 7-ray spectrum from 
Eq. 6. This condition is most directly fulfilled by the power 
law-spectrum with 7 = 2 that has equal DEBRA power 
per unit logarithmic bandwidth in energy. It results in 
t(E) oc E and then yields an upper limit for the DEBRA 
density close to e 2 n(e) = 10 _3 eV/cni 3 . This limit corre- 
sponds to the borderline on which the source spectrum 
becomes a pure power-law Figs. |l3| and [u], curves 2. A 
slight increase of the DEBRA density by as little as a fac- 
tor of 1.5 leads to a dramatic (exponential) deviation of 
the reconstructed spectrum at the highest observed 7-ray 
energies around 16 TeV from the E~ 19 power-law extrap- 
olation. 

A power law with 7 = 3 would give similar and com- 
plementary results. The resulting source spectra (Fig. 14) 
and the upper limits on the DEBRA density (Fig. 13) 
obtained in this way assuming power-law spectra for DE- 
BRA with 7 = 2 and 3, are given by the curves 2 and 
3, respectively. The power law 7=1 corresponding to 
t 77 = 3 complements Fig. |l3[ 

It should be noted however that any realistically ex- 
pected spectrum of the DEBRA in a broad range of wave- 
lengths deviates from a simple power-law. In fact, all 
models of the DEBRA, independently of the details, pre- 
dict two pronounced peaks in the spectrum at 1 /jm and 
100 /im contributed by the emission of the stars and of the 
interstellar dust, respectively, and a relatively flat 'valley' 
at mid IR wavelengths around 10 /im (see e.g. Dwek et 
al. 1998 ). The strong impact of the DEBRA spectrum on 



100 



calculations of the opacity of the intergalactic medium has 
been emphasized by Dwek & Slavin (1994) and Macminn 
and Primack (1996). 

Note that the power law with 7 ~ 1 characterizes the 
shape of the spectrum of DEBRA at near IR wavelengths, 
typically between 1 and several microns, and the power- 
law with 7 > 3 characterizes the DEBRA between 10 and 
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Fig. 15. The same as in Fig. [l3|, but for different param- 
eters of the power-law DEBRA. Curve 1: 7 = 1 with an 
absolute flux corresponding to r 77 = 1; curve 2: 7 = 1 
with e 2 n(e) = 3 • 10~ 4 eV/cm 3 ; curve 3: 7 = 4; curve 4: 
sum of spectra 1 and 3, truncated at 1 and 80 /im. 
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Fig. 16. The same as Fig. [l4]but for DEBRA fluxes shown 
in Fig. |l5|, with the addition of a reconstructed spectrum 
(curve 4) corresponding to curve 4 in Fig. EM 



This is seen in Fig. [13| and Fig. [bj where the measured 
fluxes or estimated upper and lower limits obtained di- 
rectly at different wavelengths of the DEBRA are shown. 

Finally, an interesting numerical criterion for the 
derivation of u pper limits on the DEBRA was suggested 
by Biller et al. 119981) . It relies on independent 7-ray obser- 
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than oc E within the observed energy range, is shown 
by the horizontal bars in Fig. [l3| 

The results described above could be "improved" as- 
suming a more realistic, r 77 = 1 optical depth for the 
7=1 power law branch at short wavelengths, and a 
steeper, 7 = 4, power law branch at long wavelengths 
(see Fig. [l5]). The latter choice for 7 is due to the rapid 
rise of the data points towards far infrared wavelengths 
A in order to fit the recent measurements of the flux at 
140 /mi by DIRBE aboard the COBE satellite. For the far 
infrared (FIR) branch the absolute flux of the power-law 
with 7 = 4 is chosen again from the condition that the dif- 
ferential 7-ray source spectrum should not exponentially 
rise at energies up to E ~ 16 TeV (see Fig. fi~6]). Note also 



that the criterion of t~ 



1 for the 7=1 branch at 



short wavelengths is pretty close to the level of the flux of 
the recent tenta tive detection of DEBRA at 3.5 /mi (Dwek 
& Arendt |l998|) . The sum of the NIR and FIR power-law 
branches with the above indices and absolute fluxes results 
in a deeper mid IR "valley" and predicts a very steep spec- 
trum of DEBRA between 30 and 100 /im. In Figs. 15 and 
|l6| this spectrum has been truncated at A ~ 80/mi that 
corresponds to the kinematic threshold of pair production 
at interactions with the maximum observed energy of 7- 
rays of about 20 TeV. The spectrum is also truncated at 
A = 1/im in order to avoid significant excess compared 
with the fluxes at optical/UV wavelengths recently de- 
rived from the Hubble Deep Field analysis (Pozzetti et al. 
19981 ). 



The effect of "reconstruction" of the 7-ray spectra of 
Mkn 501 corresponding to this "best estimate" of DEBRA 
between 1 and 80 /im is illustrated in Fig. |l6|. It shows, in 
the simple absorption picture (using two truncated power- 
laws) that even a conservative choice for the DEBRA field 
implies intergalactic extinction at all observed energies. 
Especially the reconstructed spectrum around 1 TeV could 
be considerably harder than the observed spectrum, with 
a maximum of E 2 dN/AE at 2 TeV (see Fig. g). This also 
demonstrates that it would be dangerous to interpret the 
observed spectral slope at low energies in terms of a power 
law extending from still lower energies. 



For a primary 7-ray spectrum dN/dE harder than 
E~ 2 , extending to energies E ^> ITeV, the cascade spec- 
trum at TeV energies could strongly dominate over the 
primary 7-ray spectrum. In addition, the spectrum of the 
cascade 7-rays that reach the observer has a standard 
form independent of the primary source spectrum with 
a characteristic photon index of 1.8 — 2.0 at energies be- 
tween ^100 GeV and an exponential cutoff determined 
by the condition r(d,E) = 1. Thus, somewhat surpris- 
ingly, the measured time- averaged spectrum of Mkn 501 
can in principle be fitted by a cascade 7-ray spectrum 
for a reasonable DEBRA flux level of about 10~ 3 eV/cm 3 , 
provided that the invisible source spectrum extends well 
beyond 25 TeV, say up to 100 TeV. 

For an intergalactic magnetic field (IGMF) B > 
10 -12 G the cascade 7-rays could be observed in the form 
of an extended emission from a giant pair halo with a 
radius up to several degrees formed around the central 
source (Aharonian et al. 1994 ). Although the possible ex- 
tinction of 7-rays from Mkn 501, at least above 10 TeV, 
unavoidably implies the formation of a pair halo, the TeV 
radiation of Mkn 501 cannot be attributed to such a halo 
simply due to arguments based on the detected angular 
size and the time variability of the radiation. 

However, the speculative assumption of an extremely 
low IGMF still allows an interpretation of the observed 
TeV 7-rays of Mkn 501 within the hypothesis of a cas- 
cade origin. Instead of extended and persistent halo ra- 
diation, we expect in this case that the cascade 7-rays 
penetrate from the source to the observer almost on a 
straight line. Yet at cosmological distances to the source 
even very small deflections of the cascade electrons by the 
IGMF lead to significant time delays of the arriving 7- 
rays: A t B ~ 2.4(d/l Gpc)(J5/l TeV)~ 2 (£?/10~ 18 G) 2 days 
(Plaga 1995). This implies that in order to see more or 
less synchronous activity (within several days or less) of 
Mkn 501 at different wavelengths, as was observed dur- 
ing multiw avelen gth observations of the source (see e.g. 

18 G. 



Pian et al. |1998|), we would have to require B < 10 
Although indeed quite speculative, such small magnetic 
fields on spatial scales large compared to 1 Mpc cannot 



be a priori excluded (e.g. Kronberg 1996) 



6.2.2. The effect of cascading in the DEBRA 



7. Conclusions 



The discussion of intergalactic 7 — 7 absorption effects 
is in principle incomplete without considering secondary 
radiations. Briefly, when a 7-ray is absorbed by pair pro- 
duction, its energy is not lost. The secondary electron- 
positron pairs create new 7-rays via inverse Compton 
scattering on the 2.7 K MBR. The new 7-rays produce 
more pairs, and t hus an electromagnetic cas cade d evelops 
(Berezinsk y et a l. 1990| , Protheroe & Stanev 1993 , Aharo- 
nian et al. Il994|). In fact, in our discussion of absorption 



In this paper we have presented the 1997 Mkn 501 time 
averaged spectrum as measured with the HEGRA IACT 
system in the energy range from 500 GeV to 24 TeV. The 
absence of strong temporal evolution as well as a signifi- 
cant correlation of the emission strength and the spectral 
shape in the energy region from 500 GeV to 15 TeV made 
the determination of a time averaged spectrum astrophys- 
ically meaningful. Due to unprecedented 7-ray statistics 
and the 20% energy resolution of the instrument, it was 
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measure a smooth, curved energy spectrum deeply into the 
exponential regime. We found, that the spectrum above 
0.5 TeV is well described by a power law with an expo- 
nential cutoff dN/dE = N E~ lm cxp (-E/6.2 TeV), the 
highest recorded photon energies being 16 TeV or more. 

The detection of TeV 7-rays from Mkn 501 leads to 
the unavoidable conclusion that the observed 7-radiation 
is produced in a relativistic jet with a Doppler factor 
Sj > 10. Actually, assuming a pure power-law production 
spectrum of 7-rays we may naturally explain the expo- 
nential cutoff in the observed spectrum spectrum by an 
internal 7-7 absorption in the jet. Because of the strong 
dependence of the optical depth on the jet's Doppler fac- 
tor, this hypothesis gives an accurate determination of the 
latter, Sj = 8.5. Remarkably, the uncertainty of this esti- 
mate, which is mainly due to the uncertainty in the value 
of the Hubble constant, H ~ 60^20 km/cMpc and in the 
measured energy of the exponential cutoff Eq = 6.2 TeV, 
docs not exceed 20%. However, since there could be other 
reasons for the steepening of the TeV spectrum, this es- 
timate can only be considered as a robust lower limit on 
<5j. In particular, the steepening of the 7-ray spectrum at 
the highest energies could be attributed to an exponential 
cutoff in the spectrum of accelerated particles, as well as 
- in the case of the inverse Compton origin of 7-rays - to 
the Klein-Nishina effect. 

In addition, a modification of the intrinsic (source) 
spectrum of TeV 7-rays takes place during their pas- 
sage through the intergalactic medium. The recent claims 
about tentative detections of the diffuse extragalactic 
background radiation by the DIRBE instrument aboard 
the COBE at near infrared (A = 3.5/im) and far infrared 
(A = 140^m) wavelengths, both at the ~ 10nW/cm 2 sr 
flux level imply a strong effect of intergalactic 7-7 extinc- 
tion on the observed Mkn 501 spectrum over the entire 
energy region measured by HEGRA. In particular, the 
shape of the highest energy part of the observed 7-ray 
spectrum combined with the DIRBE flux at 140/im re- 
quires a very steep {uF v oc A s with v ~ 4) spectrum of the 
DEBRA with a characteristic flux at mid infrared wave- 
lengths (A <~ 10 — 30^m) around 1 to 2 nW/m 2 ster. Due 
to the expected fiat spectrum of DEBRA at near infrared 
wavelengths (close to vF v oc A -1 ) the modification of the 
7-ray spectrum is less prominent at energies of a few TeV, 
although the absolute extinction could be very large. This 
does not allow us to draw definite conclusions about the 
absorption effect based on the analysis of the shape of the 
7-ray spectrum. Nevertheless, in the case of Mkn 501 this 
ambiguity can be significantly reduced by rather general 
arguments regarding the 7-ray energetics. Indeed, even 
relatively modest assumption about the optical depth of 
about r <~ 3 which corresponds to a flux of the DEBRA in 
the K-band (A = 2.2/xm) of about 30nW/m 2 ster creates 
uncomfortable conditions for any realistic models of the 



strong upper limit comparable with the DIRBE upper 
limit at this wavelength. 

To summarize, our excursion through the nonthermal 
physics of AGNs shows that 7-ray observations alone 
do not allow a unique interpretation of a spectrum like 
the one we have presented for Mkn 501, even though 
one can make a number of highly interesting inferences. 
In particular, our current poor knowledge about the 
intrinsic spectrum of Mkn 501 does not allow us to make 
definite conclusions about the effect of the intergalactic 
absorption of TeV 7-rays. And vice versa, the uncertainty 
in the density of DEBRA does not allow us to take 
into account the effect of the intergalactic absorption, 
and thus to "reconstruct" the intrinsic 7-ray spectrum. 
The knowledge of the latter is an obvious condition for 
the quantitative study of the acceleration and radiation 
processes in the source. Therefore we believe that decisive 
progress in this field could be achieved by the analysis of 
both the spectral and temporal characteristics of X-ray 
and TeV 7-ray emissions obtained during the multiwave- 
length campaigns of several X-ray selected BL Lac objects 
at different states of activity, and located at different 
distances within 1 Gpc. Due to the time variability of 
such sources this requires simultaneous observations 
together with extensive theoretical modeling. Given these 
we expect to obtain independent knowledge of the diffuse 
intergalactic radiation fields which can be compared with 
direct measurements and models of galaxy formation. 
We hope that the successful multiwavelength campaigns 
of Mkn 421 and Mkn 501 in 1998 with participation of 
several X-ray satellites and the HEGRA IACT system 
will provide highly interesting results in this area. 
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Appendix A: Systematic errors on the shape of 
the spectrum 

To illustrate the effect of systematic errors for determi- 
nation of energies and effective areas in more detail, we 
consider here the following model. For simplicity, we con- 
sider a detector consisting of a single Cherenkov telescope, 
which is located inside a uniformly illuminated Cherenkov 
light pool of area ^4 poo i- To good approximation, the light 
yield I (in photons per m 2 ) is proportional to the shower 
energy, 
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neglecting logarithmic corrections. The response of the 
telescope is now characterized by two quantities, namely 
the total intensity Q detected in the image (usually the 
size, measured e.g. in units of ADC channels, or, with a 
conversion factor, in units of 'photoelectrons'), and the 
signal V induced in the highest pixel(s), which is fed into 
the trigger circuitry. Both Q and V should be proportional 
to I, 



(3I;V = jI, 



(A.2) 



but they are influenced by rather different factors. While 
both j3 and 7 include the mirror reflectivity, the PMT 
quantum efficiency, and the PMT gain, V is to a much 
higher degree sensitive to the point spread function of the 
mirror, and to the shape of the signal generated by the 
PMT 0. The effective detection area is given by 



A = P(V)A P 



ool 



(A.3) 



where P is the trigger probability for a given value of V. 
The fact that triggering and energy determination are not 
based on identical quantities is the key origin of systematic 
errors, in particular in the threshold region. 

In the analysis of data, values ctmCj /?mc and 7mc are 
assumed for these constants, usually based on Monte Carlo 
simulations. The values may differ from the true values 
due to imperfections in the parameterization of the atmo- 
sphere (a) or of the optics and electronics of the telescope 
(f3, 7). The reconstructed energy E of a shower of true 
energy E is then, in the absence of fluctuations, 



E — ot MC 8 M( 



mc'-'mc^ ~~ a MC^MC^ a -^ — ■ (A-4) 

Based on eqs. | and | and using K(E',E) = 5{E' - f E) 
the rate of events with measured energy E is then 



r{E) = </>(E/f)P(cr,E/f)A pooi /f 



(A.5) 



In the analysis, the Monte Carlo simulated rate tmc is 
used to evaluate the effective area 



of the spectra in particular in the threshold region, where 
P(V) varies steeply with V. The scale factor / cannot be 
determined from IACT data alone; some external refer- 
ence is required. The second effect - the distortion of the 
spectra - disappears provided that cry// = cxmcImc, i-e., 
Jmc/Pmc = ll Pi assuming that the simulation correctly 
accounts for the statistical fluctuations determining the 
shape of P(V). This condition can be checked internally 
within the data set. In particular, one should compare the 
distribution in Q for a given value of V (or for a fixed 
range in V, V > Vq) in the data and in the simulation 
(see Sect. 4.3 and Fig. 5). If j//3 differs between data and 
simulation, the distribution of Q in the threshold region 
will be different. This comparison also tests the simulation 
of the shape of P(V). 

If the thresholds in Q are shifted by a factor 1 + e 
in the data relative to the simulation, this implies that 
Jmc/Pmc is off by the same factor, resulting in the flux 
error given by Eq. ^|. 
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